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We present a study of the decays B° Z)i*' + D*~, using 20.8 fb -1 of e + e~ annihilation data 
recorded with the BABAR detector. The analysis is conducted with a partial reconstruction technique, 
in which only the _Di*' + and the soft pion from the D*~ decay are reconstructed. We measure the 
branching fractions B(B° -> D+D* ) = (1.03 ± 0.14 ± 0.13 ± 0.26)% and B(B° -> D*+ D*~) = 
(1.97±0.15±0.30±0.49)%, where the first error is statistical, the second is systematic, and the third 
is the error due to the Df — > 4>ir + branching fraction uncertainty. From the B — > D* + D*~ angular 
distributions, we measure the fraction of longitudinal polarization Tl/V = (51.9±5.0±2.8)%, which 
is consistent with theoretical predictions based on factorization. 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 



*Also with Universita di Perugia, Perugia, Italy 



tAlso with Universita della Basilicata, Potenza, Italy 



4 



I. INTRODUCTION 

Precise knowledge of the branching fractions of exclu- 
sive B decay modes provides a test of the factorization 
approach [l|. Factorization neglects final state interac- 
tions between the quarks of the two final state mesons. 
The pattern of branching fractions for two-body B de- 
cays to modes such as D^n, D" p can be success- 
fully accommodated in such a model. However, it is pos- 
sible that the factorization assumption is not applicable 
to the decays B — > D^X, where the meson X contains 
a heavy quark. The current experimental uncertainties 
for B — * D* branching fractions Q do not allow us 

to perform a precise test of the factorization approach in 
this case. 

Further tests of factorization are provided by measur- 
ing the polarization in decays of B mesons to vector- 
vector final states. Within experimental errors, polariza- 
tion measurements are consistent with factorization pre- 
dictions for the final states D*p 0, Z>V(1450) Q, and 
D*D* @. 

In this paper we present measurements of the branch- 
ing fractions 1 B(B° — > Ds*' + D*~). We also report 
a measurement of the D* + polarization in the decay 
B° -> D* s +D*-, obtained from an angular analysis. 
These results provide tests of factorization with increased 
precision. 



II. THE BABAR DETECTOR AND DATA SET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II storage ring. An inte- 
grated luminosity of 20.8 fb _1 was recorded in 1999 and 
2000 at the T(AS) resonance, corresponding to about 22.7 
million produced BB pairs. 

A detailed description of the BABAR detector is pre- 
sented in Ref . Q • Only the components of the detector 
most relevant to this analysis are briefly described here. 
Charged particles are reconstructed with a five-layer, 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) with a helium-based gas mixture, 
placed in a 1.5 T solcnoidal field produced by a supercon- 
ducting magnet. The charged particle resolution is ap- 
proximately {5pt/pt) 2 = (0.0013p T ) 2 + (0.0045) 2 , where 
Pt is the transverse momentum given in GeV/c. The 
SVT, with a typical single- hit resolution of 10 /im, pro- 
vides measurement of the impact parameters of charged 
particle tracks in both the plane transverse to the beam 
direction and along the beam. Charged particle types 
are identified from the ionization energy loss (dE/dx) 



[1] Reference to a specific decay channel or state also implies the 
charge conjugate decay or state. The notation refers to 

either Z)j~ or Dg + . 



measured in the DCH and SVT, and the Cherenkov 
radiation detected in a ring imaging Cherenkov device 
(DIRC). Photons are identified by a CsI(Tl) electro- 
magnetic calorimeter (EMC) with an energy resolution 
a{E)/E = 0.023 • (S/GeV)- 1 / 4 © 0.019. 

III. METHOD OF PARTIAL 
RECONSTRUCTION 

In selecting candidates for the decays B ' — > D s D*~ 
with D*~ — > D°ir~ , no attempt is made to reconstruct 
the D a decays. Only the D s *^ + and the soft n~ from the 
D*~ decay are detected. In this way, the candidate selec- 
tion efficiency is higher by almost an order of magnitude 
than that obtained with full reconstruction of the final 
state. Given the four- momenta of the D s ' and ir~ , and 
assuming they originate from a -B — ► D s ' D*~ decay, 

the four-momentum of the B can be calculated up to an 

( *)-\- 

unknown azimuthal angle (ft around the Ds flight di- 
rection. This calculation uses the constraint of the known 
center-of-mass (CM) energy and the masses of the B° and 
D*~ mesons. Energy and momentum conservation then 
allows a determination of the four- momentum of the D°, 
whose square yields the (/(-dependent missing mass 

A/miss = ^(Ps - P D (.)+ - Ptt) 2 , (1) 

where P#, P D (*)+ and are the four-momenta of the 

B v , D s and the soft pion, respectively. In this analysis 
the missing mass is defined with an arbitrary choice for 
the angle (ft, such that the B° momentum p B makes the 
smallest possible angle with p w and p D (*)+ in the CM 
frame. 



IV. EVENT SELECTION 

For each event, we calculate the ratio of the second 
to the zeroth order Fox- Wolfram moments, using all ob- 
served charged tracks and neutral clusters. This ratio is 
required to be less than 0.35 in order to suppress contin- 
uum e + e~ — > qq events, where q = u,d, s, c. 

We reconstruct D+ mesons in the decay modes Df — > 
(ftn+, D+ -> K*°K+ and D+ -> K° S K+, with subsequent 
decays (ft -> K+K~, K*° -> K~n + and K° -» n+n-. 
These modes are selected since they offer the best com- 
bination of large branching fraction, good detection ef- 
ficiency, and high signal-to-background ratio. Charged 
tracks from the Df are required to originate from within 
±10 cm along the beam direction and ±1.5 cm in the 
transverse plane, and leave at least 12 hits in the DCH. 

Kaons are identified using dE / dx information from the 
SVT and DCH, as well as the Cherenkov angle and the 
number of photons measured with the DIRC. For each 
detector component d = {SVT, DCH, DIRC}, a likeli- 
hood (L^) is calculated given the kaon (pion) mass 
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hypothesis. A charged particle is classified as a "loose" 
kaon if it satisfies /L*, > 1 for at least one of the de- 
tector components. A "tight" kaon classification is made 
if the condition Yid^d /^d > 1 is satisfied. 

Three charged tracks consistent with originating from 
a common vertex are combined to form a Df candidate. 

In the case of the decay — ► (j)ir + , two oppositely 
charged tracks must be identified as kaons with both sat- 
isfying the loose criterion, and at least one, the tight 
criterion. No identification requirement is applied to the 
pion. The reconstructed invariant mass of the K + K~ 
candidates must be within 8MeV/c 2 of the nominal 4> 
mass : 8j. In the decay Df — > 4>n + , the <f> meson is 
polarized longitudinally, resulting in the kaons having a 
cos 2 6h distribution, where Oh is the angle between the 
K + and directions in the <p rest frame. We require 
cosOh\ > 0.3, which retains 97% of the signal while re- 
jecting about 30% of the background. 

In the reconstruction of the — » K*°K + mode, the 
K~ir + invariant mass is required to be within 65 MeV/c 2 
of the nominal K*° mass This wider window leads 
to a larger fraction of combinatorial background than in 
the Df — > (f>7T + mode. To reduce this background, we re- 
quire | cos0h\ > 0.5. In addition, substantial background 
arises from the decays D + — > K* q tt + and D + — > K°ir + , 
which tend to peak around the nominal Dj mass if the 
pion is misidcntificd as a kaon. This background is sup- 
pressed by requiring that the kaon daughter of the K*° 
satisfy the loose kaon identification criterion and the 
other kaon, the tight criterion. 

For the decay mode D+ -> K°K+, with K° s -> tt+tt", 
the 7r + 7r~ invariant mass must be within 15 MeV/c 2 of the 
nominal K® mass |8j , and the charged kaon must satisfy 
the tight criterion. To improve the purity of the sam- 
ple, we require the angle a between the K® momentum 
vector and the K® flight direction to satisfy cos a > 0.98. 

The invariant mass Mo s of -D+ candidates is required 
to be within three standard deviations {(Jd s ) of the signal 
distribution peak M D cak seen in the data. 

Df candidates satisfying these selection criteria are 
combined with photon candidates to form D* + — > 
candidates. Candidate photons are required to satisfy 
Ery > 50MeV, where E 1 is the photon energy in the 
laboratory frame, and E* > 110 MeV, where E* is the 
photon energy in the CM frame. When the photon can- 
didate is combined with any other photon candidate in 
the event, the pair must not form a good 7r° candidate, 
defined by a total CM energy E* > 200 MeV and an 
invariant mass 115 < A/ 77 < 155MeV/c 2 . 

The D*+ candidates must satisfy |AM - AMP cak | < 
2.5 a am, where AM poak is the peak of the signal AM = 
Al(Dfj) — M(Z?+) distribution observed in the data. 



The CM momentum of the D 
to be greater than 1.5 GeV/c 



(*)+ 



candidate is required 



D 



(*)+ 



candidates are combined with 7r candidates 
to form partially reconstructed B v — > D s D*~ can- 
didates. Since the transverse momentum of the pion in 



signal events is less than 210 MeV/c, these tracks are not 
required to have DCH hits. 

Due to the high combinatorial background in the AM 
distribution, more than one D* + tt~ candidate pair is 
found per event, with about a 20% probability from sig- 
nal Monte Carlo simulation. To select the best candidate 
in the event, the following \ 2 is calculated for each D* + 
candidate 

X 2 = [{Mi - MP oak )M] 2 + [{M Db M^/an,] 2 
+ [{AM — AM poak )/<7AAf] 2 , 

(2) 

where Mi is the measured invariant mass of the inter- 
mediate i = <fr, K*°, or K® candidate, depending on the 
-D+ decay mode, M peak is the corresponding peak of the 
signal Mi distribution, and 0$ is its width obtained from 
data. The candidate with the smallest value of x 2 in the 
event is retained. 



RESULTS 



The missing mass distributions of candidates for par- 
tially reconstructed B° — > d\*^ + D*~ decays are shown 
in Fig. n A clear signal peak is observed in all modes. 
We perform a binned maximum likelihood fit to these 
distributions. The fit function is the sum of a Gaussian 
distribution and a background function given by 



/s(M m ; ss ) 



Ci(M 



JlJ miss ) 



C 3 + {M - M miss ) 



c. 



(3) 



where Ci are parameters determined by the fit, and 
M = Mo* - Ma = 1.871 GeV/c 2 is the kinematic end 
point. The fits find 3704±232 and 1493±95 events under 
the Gaussian peak in the sum of the D+ir~ and D* + ir~ 
distributions, respectively. However, due to the presence 
of cross feed and self-cross feed, discussed below, further 
analysis is needed in order to extract the signal yields 
and the branching fractions. 

We use a Monte Carlo simulation, which includes both 
BB and qq continuum events, to study the missing mass 
distributions of the different background sources. We 
consider two kinds of backgrounds: a peaking component 
that contributes predominantly at the end of the missing 
mass distribution in the signal region and a non-peaking 
component that is more uniform. The non-peaking com- 
ponent is well modeled by the background function |J3}. 
The peaking component receives contributions from re- 
lated channels due to 

• Cross Feed (CF): if the soft photon from a 
D* + — * Df-f decay is not reconstructed, B° — » 
D* + D*~ decays may lead to an enhancement un- 
der the signal peak of the Dfn" missing mass spec- 
trum. Similarly, B° — > D+D* ~ decays may lead to 
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200 



500 



CO 
CD 



100 




100 



200 



50 



(b) 




o 

1.88 1.78 1.83 

M mi5S (GeV/c 2 ) 



TABLE I: Efficiencies for B° -» L>i* )+ D*~ decay modes to 
contribute to the Dfn~ and Dg + 7r~ missing mass distribu- 
tions in the signal region M miaa > 1.86GeV/c 2 . Two differ- 
ent B — > D* a + D*~ Monte Carlo samples have been used, 
one with longitudinal (long.) and the other with transverse 
(transv.) polarization. 





Reconstructed mode 


Generated mode 


D+7V- 




B u — > DfD*~ 


(23.6 ± 1.0)% 


(1.7 ±0.3)% 


B u -> D"+D- (long.) 
Self-Cross Feed 


(9.0 ± 0.3)% 


(7.4 ± 0.3)% 
(1.6 ±0.1)% 


B u -> D* s + D*~ (transv.) 
Self-Cross Feed 


(10.4 ±0.3)% 


(6.9 ±0.3)% 
(1.4 ±0.1)% 




FIG. 1: Missing mass distributions of B candidates in data, 
(a) D+tt~ with Df -> c/>7r_ + , (b) D* s + tv' with D+ -> 4>ir + , 
(c) _D+7r~ with D+ -» (d) DJ+tT with D+ -> 

K*°K+, (e) D^tv~ with £>+ -> A"sA: + , (f) D^+tt" with 
D+ -> A'sAT + . The curves show the result of the fit (see 
text), indicating the signal and background contributions. 



a peaking enhancement in the D* + n~ M m - lss spec- 
trum, due to the combination of a Dj with a ran- 
dom photon. 

Self-Cross Feed (SCF): this is due to true B° -> 
D* + D*~ decays in which the Df is correctly re- 
constructed, but combined with a random photon 
to produce the wrong D* + candidate, resulting in 
a peaking enhancement in the D* + ?t~ spectrum. 



Table [I] presents the reconstruction efficiency of correctly 
reconstructed signal B° — > Dg + D*~ decays, as well as 
cross feed and self-cross feed, found for simulated events 
in the signal region M miss > 1.86 GeV/c 2 . 

In addition to these background sources, we also con- 
sidered a possible contribution from the charged and neu- 
(*)+ — 

tral B decays B — > D s D** . These potential back- 
ground sources were simulated with four D** states: 
the observed -Di(2420) and £>2(2460) mesons, and the 
Do(J = V2) an d Di(j = 1/2) mesons predicted by 
HQET 9]. Their contribution was determined to be 

(*)+ 

negligible, mainly due to the D s CM momentum cut. 




1.85 1.875 
^WGeV/c 2 ) 



FIG. 2: Missing mass distribution for (a) Dfir~ and (b) 
D* + 7T~ combinations for data (points with error bars) and 
Monte Carlo (histogram). The contributions from the BB, 
cc and qq with q = u, d, s (labeled uds in the figure) are 
shown separately. The cross feed and self-cross feed back- 
grounds are included in the total histogram, but not in the 
hatched BB histogram. 



Multi-body decays B — > D 



^ + X are found not to con- 
tribute due to the same cut. 

Figure |21 shows a comparison of the missing mass dis- 
tributions in data and Monte Carlo simulation. We 
assume 1.05% and 1.59% branching fractions for the 
B° -> D+D*- and B° -> D*+D*~ decays, respectively, 
in the Monte Carlo simulation. 
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The number of events in the peaks in the Dfn~ and 
D* + 7r - Af m i SS distributions is obtained from the fits de- 
scribed above. The branching fractions are computed 
from these yields correcting for cross feed and self-cross 
feed background. This is done by inverting the 2 x 2 ef- 
ficiency matrix, whose diagonal elements correspond to 
the sum of signal and self-cross feed efficiencies presented 
in Table [I] and whose off-diagonal terms are the cross- 
feed efficiencies. The efficiencies corresponding to trans- 
verse and longitudinal polarization of B° — > D* + D*~ 
have been weighted according to the measured polariza- 
tion discussed below. With this procedure, the B° — > 
branching fractions are determined to be 

D+D*-) = (1.03 ± 0.14 ± 0.13 ± 0.26)%, (4) 



D\* )+ D 
B{B° 



B(B° -> D* + D*-) = (1.97 ± 0.15 ± 0.30 ± 0.49)%, (5) 

and their sum is 

B{B° -> D^+D*-) = (3.00 ±0.19 ±0.39 ±0.75)%, (6) 

where the first error is statistical, the second is the sys- 
tematic error from all sources other than the uncertainty 



in the D: 



i>7r + branching fraction, and the third 



error, which is dominant, is due the uncertainty in the 
Df — > 4>tt + branching fraction B(Df — > (f>ir + ) — (3.6 ± 
0.9)% 8J. Correlations in the systematic errors between 
Eqs. (@J and J5J are taken into account in Eq. ©. The 
sources of the systematic error are discussed in Sec. IVII 

The measurement of the fraction of the longitudinal 
polarization T L /T in the B a — » D* + D*~ decay mode 
is performed for candidates having missing mass in the 
signal region (M m i ss > 1.86GeV/c 2 ). To minimize the 
systematic error due to large backgrounds, the polariza- 
tion measurement involves only the — > 07r + channel, 
which has the best signal to background ratio. Two an- 
gles are used: the helicity angle 6 7 between the D*~ and 
the soft photon direction in the D* + rest frame, and the 
helicity angle 6^ between the D* + and the soft pion di- 
rection in the D*~ rest frame. Since the B meson is not 
fully reconstructed, we compute 6 1 and 0^ by constrain- 
ing Mmiss to the nominal D° mass 8] to obtain a unique 
solution for the azimuth <f>. 

The two-dimensional distribution (cos# 7 , cos 9^) is di- 
vided into five ranges in each dimension, resulting in 25 
bins. The combinatorial background, as well as the cross 
feed and the self-cross feed obtained from Monte Carlo 
simulation, are subtracted from this two-dimensional 
data distribution. The resulting signal distribution is 
corrected bin-by-bin for detection efficiency, which is 
obtained from the simulation separately for each bin. 
A two-dimensional binned minimum-x 2 fit is then per- 
formed to the efficiency-corrected signal distribution with 
the fit function 



d 2 Y 



dcos^ d cos ft, 



r 

(i - 



cos 



Sill 




FIG. 3: Projections of the number of background-subtracted 
data events on the cos6>,r and cos# 7 axes. The result of the 
two-dimensional fit is overlaid. 



The resulting fit has a % 2 of 23.1 for 25 bins with two 
floating parameters (r^, /Y and total normalization) . Fig- 
ure |3| shows the data and the result of the fit projected 
on the cos# 7 and cos#,r axes. 

From the fit, the fraction of longitudinal polarization 
is determined to be 



Y L /Y = (51.9 ±5.0 ±2.8)%, 



(8) 



Y 



'-(7) 



where the first error is statistical and the second is sys- 
tematic. 



VI. SYSTEMATIC UNCERTAINTIES 

The various contributions to the systematic errors on 
the branching fraction and polarization measurements 
are summarized in Table [H] The dominant systematic 
error for the branching fractions is the uncertainty on 
the three Dj branching fractions. To evaluate the uncer- 
tainty due to the background subtraction, the signal yield 
is determined using an alternative method in which the 
number of events is extracted directly from the histogram 
after subtraction of the background, which is estimated 
with the Monte Carlo simulation. The difference of the 
signal yields obtained in this way from the results of the 
fit was taken as a systematic error. This also accounts 
for the systematic error due to a possible deviation of the 
signal shape from a Gaussian. 

The Monte Carlo statistical errors in the determination 
of the signal and the cross feed efficiencies are included 
as a systematic error. The uncertainty in the calculation 
of the B° — > D* + D*~ polarization is propagated to the 
branching fraction systematic error. The systematic er- 
ror due to the uncertainty on the efficiency for the recon- 
struction of charged particles is 1.2% times the number 
of charged particles in the decay. An additional error of 
1.6% is added in quadrature to account for the uncer- 
tainty in the reconstruction efficiency of the soft pion. 



8 



TABLE II: Sources of systematic uncertainties (%) for the 
B° -> Z)i* )+ D*- branching fractions and B° -> D* + D*~ 
polarization measurements. 



Source 


D+D*~ 


D*+D*~ 


*(Ti/r) 


Background subtraction 








or modeling 


2.7 


5.9 


0.5 


Monte Carlo statistics 


4.2 


6.0 


2.7 


Polarization uncertainty 


0.8 


0.5 


- 


Cross feed 


3.2 


2.4 


- 


Number of B pairs 


1.6 


1.6 




B(<f> -» K+K-) 


1.6 


1.6 




Particle identification 


1.0 


1.0 


0.1 


Tracking efficiency 


3.6 


3.6 


0.5 


Soft pion efficiency 


2.0 


2.0 


0.2 


Relative branching fractions 


10.2 


10.2 




B(D* S + -> D+ 7 ) 




2.7 




Photon efficiency 




1.3 


0.1 


7T° veto 




2.7 


0.3 


Total systematic error 


13.1 


15.1 


2.8 



The systematic error due to the ir° veto requirement was 
studied by measuring the relative D* + yields in data and 
Monte Carlo with and without the ir° veto. 

In the polarization measurement, the level of the var- 
ious backgrounds depends on the charged track, neutral 
cluster, and particle identification efficiencies. The fit 
was repeated varying the background according to the 
errors in these efficiencies, and the resulting variations in 
Tl /r were taken as the associated systematic error. 



TABLE III: The fractional difference ((JV D - N M c)/N M c), 
averaged over all M m i ss bins, where Ad (Nmc) is the num- 
ber of data (Monte Carlo) candidates in a given bin of the 
Mmiss distribution of the given control sample. SB (SR) refers 
to the Md s or AM sideband (signal region) control sam- 
ple. WS indicates wrong sign Dl* ' + 7r + combinations, and 
—p* D („)+ indicates that A/ m i ss was calculated from the nega- 
tive of the Di*^ + CM momentum. The missing mass range 
1.78 < A/ m i ss < 1.87GeV/c 2 is used for the control sample, 
except for the first line. 



Sample type 



Dtn~ 



1.78 < M, 
SB 

SR, WS 
SB, WS 
SR, -p" 
SB, - 



< 1.85GeV/c 2 



-0.009 ± 0.007 
-0.075 ± 0.006 

0.006 ± 0.008 
-0.060 ± 0.007 

0.015 ±0.009 



0.075 ±0.014 
0.007 ±0.022 
0.044 ±0.015 
-0.008 ± 0.024 
0.075 ± 0.016 



-0.062 ±0.007 -0.123 ±0.022 



Average 



-0.038 ± 0.003 0.032 ± 0.007 



To check that the simulation accurately reproduces 



the background M m ; ss distributions in the data, a thor- 
ough data-Monte Carlo comparison is made in control 
samples containing no signal events. These samples 
are events with 1.78 < M m i SS < 1.85GeV/c 2 ; events 
in the D+ sideband 1.89 < M D , < 1.95GeV/c 2 or 
1.985 < M Ds < 2.05GeV/c 2 ; events in the D* s + sideband 

170 < AM < 300MeV/c 2 ; wrong sign Z)i* )+ 7r + combina- 
tions in either the Md s and AM sidebands or signal re- 
gions determined above; and candidates in which M m ; ss 
was calculated with the inverse of the D y s ' center-of- 
mass momentum p 



The comparison between data 
and Monte Carlo simulation for these control samples is 
shown in Table II I II The average level of discrepancy is 
used to estimate the uncertainty in the modeling of the 
background. 



VII. SUMMARY 

In summary, based on a partial reconstruction tech- 
nique, we have measured the branching fractions 

B(B° D+D*-) = (1.03 ± 0.14 ± 0.13 ± 0.26)% , 
B(B° D*+D*-) = (1.97 ± 0.15 ± 0.30 ± 0.49)% , 
B{B° -> Di* )+ D*-) = (3.00 ± 0.19 ± 0.39 ± 0.75)% . 

The fraction of the longitudinal D* + polarization in 
B° —> D* + D*~ decays is determined to be 

T L /r= (51.9 ±5.0 ±2.8)%. 

This measurement is consistent with theoretical predic- 
tions assuming factorization, which range from 50 to 
55% 03 E3. Our results are also in good agreement 
with previous experimental results 
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